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Abstract 

We present a complete analysis of the cosmological constraints on decaying dark matter. Previous 
analyses have used the cosmic microwave background and Type la supernova. We have updated 
them with the latest data as well as extended the analysis with the inclusion of Lyman-a forest, 
large scale structure and weak lensing observations. Astrophysical constraints are not considered 
in the present paper. The bounds on the lifetime of decaying dark matter are dominated by 
either the late-time integrated Sachs- Wolfe effect for the scenario with weak reionization, or CMB 
polarization observations when there is significant reionization. For the respective scenarios, the 
lifetimes for decaying dark matter are P"^ > 100 Gyr and (/P)~^ > 5.3 x 10^ Gyr (at 95.4% 
confidence level), where the phenomenological parameter / is the fraction of the decay energy 
deposited in baryonic gas. This allows us to constrain particle physics models with dark matter 
candidates through investigation of dark matter decays into Standard Model particles via effective 
operators. For decaying dark matter of ~ 100 GeV mass, we found that the size of the coupling 
constant in the effective dimension-4 operators responsible for dark matter decay has to generically 
be < 10~^^. We have also explored the implications of our analysis for representative models in 
theories of gauge-mediated supersymmetry breaking, minimal supergravity and little Higgs. 
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1. INTRODUCTION 



The past decade and a half has seen tremendous progress in the field of cosmology. The 
numerous experiments and the plethora of accumulated data have elevated cosmology into 
a precision science. With the knowledge we have gained, a remarkably consistent consensus 
known as the standard model of cosmology has emerged from these attempts to understand 
the nature of the universe. The picture that we have presently is that of a universe that 
is composed of 74% dark energy, 22% dark matter (DM) and 4% baryonic matter [H, 
Despite its extraordinary success in explaining a variety of diverse observations, fundamental 
questions do remain. Arguably the most vexatious is the question, "What are these dark 
components of the universe?" Despite the fact that it makes up 96% of the universe, we have 
so far been unable to say definitively what they are. Perhaps the most compelling ideas to 
resolve this question arise from particle physics. Dark energy is commonly attributed to 
the vacuum energy while the identity of dark matter is hypothesized to be one of the new 
particles in theories that extend the Standard Model of particle physics. This confiuence 
and cross-fertilization of ideas from two major fields of scientific endeavor promises to herald 
an exciting new era in understanding of the universe. With the resumption of operation of 
the Large Hadron Collider at CERN, we are possibly months away from collider detection, 
albeit indirectly through missing energy signatures, of dark matter particles. 

Yet this avenue too is plagued with troubling questions that need to be addressed if we 
are to take the idea seriously that some particle from an extension of the Standard Model 
is indeed the elusive yet ubiquitous dark matter of the universe. Most worryingly, to ensure 
the presence of a dark matter candidate in a number of beyond the Standard Model theories 
of particle physics, it is often necessary to impose global symmetries. For instance, we have 
the T-parity in Little Higgs and R-parity in supersymmetry 0]. In the limit where the 
global symmetries are exact, the lightest particle carrying such a global charge would be 
stable from decay to lighter particles that do not possess such a charge. It is this point 
that could potentially destroy this promising marriage of ideas from cosmology and particle 
physics, for it is well known that global symmetries are never exact. 

The presence of anomalies, as in the case of T-parity [sj], or R-parity violating terms 
in supersymmetry would often mean that the dark matter candidates arising from these 
theories are neither stable nor long-lived in the cosmological sense. Even if this had not been 
the case, the presence of gravity would necessarily induce the violation of global symmetries 
as was first revealed in studies of black holes 0, S, 0, [13]. So the lightest particle charged 
under a particular global symmetry would have, at best, a very long lifetime. Indeed, it has 
even been conjectured that discrete global symmetries are violated maximally by gravity 
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Additional motivation can be found in numerical simulations of the universe (based on 
the conventional ACDM cosmology) which predict an overabundance of substructures as 
compared to actual observations. Models with decaying dark matter [l^, l3] provide an 



extremely compelling and natural mechanism for suppressing the power spectrum at small 
scales thus resolving the discrepancy. 

Continuing on the line of thought leading from particle physics to cosmology, the question 
that then naturally springs to mind is 'What can cosmology say about decaying dark matter 
and the particle physics theories that contain them?" It is this intriguing prospect that we 
explore in this paper. 



There have been a few papers [15|, [la, [iTj in recent years analyzing DM decay into 
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electromagnetically non-interacting particles using just the cosmic microwave background 
data (Ref.[l6[] also includes supernova data). In this paper, we revisit this scenario using 
a Markov Chain Monte Carlo (MCMC) analysis employing all available datasets from the 
cosmic microwave background (CMB), Type la supernova (SN), Lyman-a forest (Lya), large 
scale structure (LSS) and weak lensing (WL) observations. We find that the hfetime of 
decaying DM is constrained predominantly by the late time Integrated Sachs Wolfe (ISW) 
effect to be F"^ > 100 Gyr. In the main body of the paper, we will comment on the 



discrepancies between the results of Refs. [15|, lla, III 



The studies in the preceding paragraph considered only the case where there was negli 
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gible reionization of the universe due to DM decay. In an attempt to address this, Ref. 
analyzed the scenario of DM decaying into only electromagnetically interacting products, 
that get partially absorbed by the baryonic gas, using a subset of the available CMB datasets. 
Our paper extends their analysis by using all the available CMB datasets, and also the SN, 
Lya, LSS and WL datasets. Besides the smaller selection of datasets, their analysis also 
ignores the impact of DM decay on cosmological perturbations which renders it ineffectual 
in the parameter space where there is negligible reionization. Our treatment allows the 
decay products to not only reionize the universe but also takes into account the effect of 
DM decay on cosmological perturbation. This allows us to generate many other observ- 
ables, particularly, the late time ISW effect that is crucial to constrain the lifetimes at low 
reionization. Another key difference between the analyses is that we use a combined reioniza- 
tion parameter for both DM reionization and phenomenological star formation reionization, 
rather than just treating them separately as was done in Ref.[l8|, because current observa- 
tions cannot distinguish which contribution to reionization is the dominant one. Doing the 
MCMC analysis, we find that the lifetime of decaying DM in this scenario constrained to be 
(/r)~^ > 5.3 X 10* Gyr, where / is a phenomenological parameter introduced by Ref. 1^ 
and related to the degree of reionization. 

Astrophysical constraints together with additional assumptions have also been used 



|19l . |20| to give even tighter bounds than ours on the lifetime of the decaying DM. While 
interesting and complementary, these lie outside the scope of our present paper. 

Having obtained the bounds on the lifetime of decaying dark matter, we will then explore 
the implications of our cosmological analysis on particle physics models beyond the Standard 
Model. We will present a complete list of cross-sections for spin-0, spin-1/2 and spin-1 dark 
matter to decay into Standard Model degrees of freedom via effective operators. Obviously, 
this can be easily extended to other models with additional light degrees of freedom (for 
instance, hidden valley models [2l|) by appropriate substitution of the parameters. Applying 
the bound on the lifetime of the decaying DM, we can then place limits on the size of the 
parameters of theories. For generic theories with a decaying dark matter of ~ 100 GeV mass, 
the coupling constant in the effective dimension-4 operators responsible for dark matter 
decay will be shown to be < 10~^^. We will also look at specific representative cases of 
theories beyond the Standard Model physics and investigate the possibility of viable dark 
matter candidates: the spin-0 messenger DM in the context of gauge mediation messenger 
number violation, the spin-1/2 bino DM in the scenario with R-parity violation and the 
spin-1 "massive photon partner" DM in the framework of T-parity violation. 

The rest of this paper is organized as follows. In Section 2, we discuss the physics 
of decaying dark matter cosmology as well as introduce the datasets that we will be using. 
Section 3 contains our Markov Chain Monte Carlo results and discussions of the cosmological 
implications. In Section 4, we explore the consequences of these results for particle physics 
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theories by enumerating the decay channels and partial widths. Representative models from 
theories of gauge-mediated supersymmetry breaking, minimal supergravity and little Higgs 
were also investigated using the results of our analysis. We conclude and briefly comment 
on future prospects in Section 5. 



2. DECAYING COLD DARK MATTER COSMOLOGY 

We will assume the standard picture of ACDM cosmology, i.e. a Friedman-Robertson- 
Walker universe that is principally composed of dark energy and cold dark matter, with 
one crucial modification; that is, we have a cold dark matter that is very long-lived but 
ultimately decays. As we are considering lifetimes of gigayears (Gyr), the fraction of DM 
decays happening during or soon after big bang nucleosynthesis (BBN) is negligible and 
hence would not alter the predictions of BBN. To perform a model-independent analysis, 
we allowed decays to all possible SM particles. However, we will assume that the long 
term decay products are relativistic. While we include branching ratios to intermediate 
non-relativistic states, they are assumed to be short-lived and will rapidly decay into light 
relativistic degrees of freedom. 

The evolution of background and first order perturbation in decaying cold dark matter 
model was first formulated in longitudinal gauge [l5|, which means the decay rate has to 
be treated with care as the CDM is not at rest in the longitudinal gauge. We, on the other 
hand, will work in CDM rest frame using synchronous gauge with the line element written 
as 

ds^ = a^(r)[-dr^ + {5ij + hij)dx' dx^]. (1) 

where r is conformal time, and t the cosmological time {dt = a{T)dr). In this paper we 
follow the convention a = 1 today. 
The decay equation 

dpcdm 

dt 

can be reformulated in a covariant form 



-^Pcdrrn (2) 



n;^(CDM) = G„ (3) 
where the force density vector Giy can be calculated from its value in CDM rest frame 

G'j^lcDM rest — 

(-^Ped„^, 0,0,0). (4) 

The conservation of total energy momentum tensor requires 

n^dr) = -G., (5) 

where the daughter radiation (dr) is composed of the CDM decay products. 
The equations describing the evolution of background are 

pcdm = —'iT^Pcdm — Pcdm, (6) 
Pdr = -4:T-Cpdr + a^Pcdm, (7) 

where dot denotes the derivative with respect to conformal time r. We have defined the 
conformal expansion rate to be 7i = ^. 
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We will only consider scalar metric perturbations, which in Fourier space can be expanded 
as following 2^ • 



/iii(x, = j rf^ke*'''''[ninj/i(k, r) + 6(n,nj - l6ij)r]{k, r)], (8) 
where n = k/|k|. 

Our choice of gauge and coordinates lead to the following simple density perturbation 
equation for CDM, 

The terms containing F all cancel out, because the background density and overdensity are 
decaying with the same rate. 

Instead of using simple hydrodynamic approximation for the decay product jlsj, which 
might give correct order of magnitude but less accurate results, we use the full Boltzmann 
equations to describe the decay product, which were first given by Ref.jl^], and recently 
updated by [12] for decaying DM cosmology. 



Ih- ^kvdr + aT(^^{6cdm-Sdr), (10) 

Pdr 



'Jdr — 

Vdr = k{\6dr - lUdr) " aT^Vdr, (H) 

tldr = Kl^Vdr - If,) + ±h + li^- aF^n,,, (12) 

Pdr 



'I 



^ -[/F,_i-(/ + l)F,+i]-aF^F,, (13) 



2/ + 1 pdr 



where / = 3,4,5,...,, F2 = 11 and for the rest, we have used the conventions of Ref.[22[. 
Because CDM particles are heavy and non-relativistic, we have treated the CDM as a perfect 
fluid. 

For the case where the DM candidate also decays into electromagnetically interacting 
particles (e.g. photons or electron/positrons), we have to be more careful. This is because 
the decays may deposit significant energy into baryonic gas and contribute to the reionization 



of universe. Following [18|, |25|, |26[ , we introduced a phenomenological factor / as the fraction 
of the decay energy deposited in the baryonic gas. For long-lifetime dark matter models, 
the reionization due to dark matter decay only depends on the combination C, = fT/Ho. 
We use ( as an additional parameter in our MCMC analysis. Without a prior on /, the 
constraint on ( does not directly give any information on F. However, for given dark matter 
models, one should in principle be able to calculate the decay branches, and therefore give 
a rough estimate for f under certain additional assumptions. Following Ref. 18, 2^, we 
modify RECFAST 27, 2^ to calculate the reionization due to DM decay. In this scenario, 
the reionization is dominated by DM decay at redshift z > 20, and is competing with the 
contribution from star formation (or other sources) at some redshift between z = 6 and 
z = 20. Without knowing the details of star formation or other reionization sources, we use 
the following phenomenological model, which can be regarded as a combination of CosmoMC 
phenomenological formula and DM decay reionization formula, 

= maxjxe^^^^^^^, ^-^[l + tanh(^i±^^^i^^ii^^)]}, (14) 
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where Xe, the ionized fraction, is defined as the ratio of free electron number density to hydro- 
gen number density; fne is the ratio of hehum number density to hydrogen number density; 



X. 



RECFAST 



is the modified RECFAST output ionized fraction (i.e., the ionized fraction as- 
suming DM decay is the only source of reionization; Az is the redshift width of reionization 
(due to other sources), for which we have taken the CosmoMC j29| default value 0.5; the 
last free parameter, reionization redshift Zre, is determined by the total optical depth Tre- 

With all the above equations, we modified CosmoMC to analyze the decaying CDM 
model. In addition, we also incorporated weak lensing data into the Markov Chain Monte 
Carlo (MCMC) analysis. The datasets used in this paper are listed below. For each dataset, 
we either wrote a new module to calculate the likelihood or modified the default CosmoMC 
likelihood codes to include the features of the decaying CDM model. 

Cosmic Microwave Background ( CMB ) 



We employ the CMB datasets from WMAP-5yr 



BOOMERANG 30, 31, 32 



ACBAR 
li 



34, 35, 36 



CBI HHHIij, VSA [iil, DASI Hi!, and MAXIMA 
Also included are the Sunyaev-Zeldovich (SZ) effect for WMAP-5yr, ACBAR 
and CBI datasets. The SZ template is obtained from hydrodynamical simulation 45 



When calculating the theoretical CMB power spectrum, we have also turned on CMB 
lensing in CosmoMC. 

Type la Supernova (SN) 

The Union Supernova la data (307 SN la samples) from The Supernova Cosmology 
Project [IgI was utilized. For parameter estimation, systematic errors were always 
included. 



Large Scale Structure (LSS) 



For large scale structure we will use the combination of 2dFGRS dataset [47|] and SDSS 



Luminous Red Galaxy Samples from SDSS data release 4 [48] . It should be noted that 



the power spectrum likelihood already contains the information about BAO (Baryon 



Acoustic Oscillation [491. |50[|). 



Weak Lensing (WL) 

Five weak lensing datasets were employed in this paper. The effective survey area and 
galaxy number density of each survey are listed in Table [H 



TABLE I: Weak Lensing Data 



Data 


Ae//(deg-2) 


ne//(arcmin ^) 


COSMOS [51, 52] 


1.6 


40 


CFHTLS-wide [53, 54, 55] 


22 


12 


GaBODS [54, 56, 57] 


13 


12.5 


RCS [54, 58, 59] 


53 


8 


VIRMOS-DESCART [54, 55, 60] 


8.5 


15 



For COSMOS data we used the CosmoMC module written by Julien Lesgourgues|52l| . 
For the other four weak lensing datasets, we utilized the likelihood given in Benjamin 
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et al. [H^]. We take the best fit parameters for tlie following galaxy number density 
formula, 



n[z] 



-rexp (--)''. (15) 



Benjamin et al. give two sets of best fit parameters, fitting on the galaxy samples 
with median photometric redshift < < 4 and 0.2 < Zp < 1.5, respectively (see 
Fig. 2 and Table 2 in [31 )• We only used the data from the Q < Zp < A region. We 
also simplified the marginalization on n{z) parameters by assuming a Gaussian prior 
on zq. The width of Gaussian prior is adjusted so that the mean redshift Zm has an 
uncertainty of 0.03(1 + 2^), i.e. 

r(i±^) 

^.o = 0.03[^;^ + ^o]. (16) 

The weak lensing data only measure matter power spectrum at angular scales less than 
a few degrees, which corresponds to scales less than a few hundred Mpc. This is much 
less than the Jean's length of the daughter radiation and therefore we can ignore the 
daughter radiation when calculating the power spectrum of projected density field. 

We should stress that this is specific to the decaying CDM model and differs from the 



conventional CDM model l61l. 62 



Lyman-a Forest 

The following Lya forest datasets were applied. 

1. The dataset from Viel et al. ^] consist of LUQAS sample ||64] and the Croft 
et al. data [65l. 



2. The SDSS Lya data presented in McDonald et al. [g^, 6^. To calculate the 
likelihood, we interpolated the table in the three-dimensional amplitude-index- 
running space. 

We explore the likelihood in nine- dimensional parameter space, i.e., the Hubble parameter 
h, the baryon density flbh'^, the amplitude and index of primordial power spectrum ( Ag and 
Tig), the DM decay reionization parameter (, the total reionization optical depth r^e, the SZ 
amplitude Asz, the decay rate normalized by Hubble parameter and the CDM density 
in early universe flcdm,eh'^- The parameter flcdm,e is defined to be 

^ _ (PcdmQ^)|a<l 

PcritO 

where pcrito = |^ is today's critical density. As the CDM in our case decays, we made a 
distinction between flcdm,e and ^cdm where the latter is defined to be the usual fractional 
CDM density today {pcdmo / Pcrito) ■ 
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FIG. 1: Posterior probability density function of the decay rate T. Solid line: using all the datasets. 
Dashed line: CMB + SN + LSS + Lya. Dotted line: CMB only. The probability density function 
is normalized as J P{T)dT = 1. 




0.005 0.01 0.015 0.005 0.01 0.015 

r/Gyr ' r/Gyr ' 



FIG. 2: Constraints on the early universe CDM density parameter i^cdm,e and decay rate F, using all 
the datasets, is plotted on the left panel. For comparison, the present day CDM density parameter 
^cdm and decay rate F is plotted on the right panel. The inner and outer contours correspond to 
68.3% and 95.4% confidence levels, respectively. 



3. MARKOV CHAIN MONTE CARLO RESULTS AND DISCUSSION 

For the case with negligible reionization, we generated 8 MCMC chains, each of which 
contains about 3000 samples. The posterior probability density function of CDM decay rate 
can be directly calculated from the Markov Chains, as shown in Fig. [T] The corresponding 
68.3% and 95.4% confidence level lower bounds on lifetime are F"^ > 230Gyr and F~^ > 
lOOGyr, respectively. If we take the lifetime of universe to be 14Gyr, the 95.4% confidence 
level limit (i.e. lifetime lOOGyr) corresponds to a scenario that roughly 15% of CDM has 
decayed into radiation by today. 

In our analysis, all the early universe physics before recombination remains unchanged. 



8 




FIG. 3: CMB power spectrum for different dark matter decay rate, assuming tlie decayed particles 
are relativistic and weakly interacting. For the CDM density parameter, we choose ^cdm,eh'^ to be 
the same as WMAP-5yr median ilm^^- For the other cosmological parameters we use WMAP-5yr 
median values. By doing this, we have fixed the CDM to baryon ratio at recombination. In a 
similar plot in Ichiki et al. 15| ilcdmoh"^ is instead fixed. Therefore the height of first peak, which 
has dependence on CDM to baryon ratio at recombination, will significantly change as one varies 
the decay rate. In this plot the red line corresponds to a stable dark matter . The blue dotted line 
corresponds to dark matter with a lifetime 100 Gyr, and the blue dashed line 27 Gyr. The data 
points are WMAP-5yr < TT > spectrum mean values and errors (including instrumental errors 
and cosmic variance). 



The CMB power spectrum is however significantly modified due to two effects. One is 
that the decay of CDM modifies the evolution of background, which results in a different 
distance to last scattering surface compared to the conventional case. The second one is that 
the decay of CDM affects the cosmological perturbations in late universe, resulting in an 
enhancement of the integrated Sachs- Wolfe (ISW) effect beyond that due to the cosmological 
constant. And it is this effect, anticipated by Kofman et al. jl^], that gives us the most 
restrictive bound on the lifetime of decaying dark matter for the scenario with negligible 
reionization. 

Let us review the inconsistency between past papers on this issue. We start with Ref . [l6| . 
Now, CMB and SN observations today can measure the fractional CDM density to a roughly 
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15% level (within 95% confidence level). We will expect the constraints on CDM decay 
ratio to be the same order of magnitude. This simple estimation does not take into account 
the fact that the decayed product still forms part of the matter component (the equation of 
state is changed to |), and that the DM decay happens mostly at low redshift. Therefore 
if we do not take the cosmological perturbation into account, the data should allow about 
15% dark matter to have decayed by today, i.e., we should not get a bound better than 100 
Gyr. This simple analysis implies that the recent lower bound of lifetime (r~^ > 700Gyr 
at 95.4% confidence level) obtained by Ref.[l6|, which does not take into consideration the 
impact of DM decay on cosmological perturbation (the location of first CMB peak is affected 
only through the change of background evolution), may not be credible. If indeed the CDM 
lifetime is 700Gyr, only about 2% of CDM has decayed into radiation by today, and by 
the time of recombination, less than 10~^ of CDM has decayed. The change in background 
evolution is so tiny that it should not be detectable by current cosmological data. 



To compare with Ichiki et al. [15|, we re-did the analysis using just the CMB datasets, 
and found the CDM lifetime > 70 Gyr at 95.4% CL, which is consistent with their 
results. The reason we obtained a more constrained value than their > 52 Gyr at 95.4% 
C.L. is probably because we used WMAP-5yr compared with their WMAP-lyr dataset. We 
expect the WMAP-9yr dataset, when published and analyzed, to exhibit only a modest 
improvement because the information from the late time ISW effect is limited by cosmic 
variance. Recently, Lattanzi et al. 17| obtained a bound of > 250 Gyr at 95.4% C.L. 
with just the WMAP-3yr data, which is not consistent with both Ichiki et al. and our 
results. We notice that in the Fig. 3 of their paper, the proximity between 68.3% and 95.4% 
confidence level bounds on F indicates a sudden drop of marginalized likelihood C{T). In 
our result, as shown in Fig. [T], this sudden drop feature is not seen. 

Let us move on to the scenario where there is significant reionization due to the decaying 
dark matter. We generated another 8 MCMC chains, each of which contains about 6000 
samples. The results of our analysis can be seen in Fig. H] and O where we show the 
constraints on DM decay reionization parameter. A few things should be noted. Firstly, the 
sharp boundary (reflected in the closeness of the two contours) on the rising edge in Fig. [5] is 
due to the fact that for a given r^e, DM decay has an upper limit because the optical depth 
due to DM decay should not extend beyond r^e. Secondly, the plateau of likelihood around 
/ F = in Fig. H] indicates that current CMB polarization data can only constrain the total 
optical depth, but cannot distinguish between DM decay reionization and star formation 
reionization. In other words, the data does not favor or disfavor DM decay reionization, as 
long as its contribution to total optical depth is not larger than the preferred r^e. 

The constraint we have obtained is / F < 0.59 x 10~^^s~^ at 95.4% confidence level. This 



result is about a factor of 3 better than Zhang et al. [18|. In the limit where reionization 
is negligible, Zhang et al. cannot give a strong bound on F because they have ignored the 
impact of DM decay on cosmological perturbations. Hence, their constraint on DM decay is 
essentially, only from CMB polarization data. Our analysis, which combines many different 
cosmological datasets and includes the calculation of the impact of DM decay on all the 
observables, gives a stringent constraint on F even in the / = limit. As for the limit of 
significant reionization, our bounds are, as mentioned earlier, an improvement over Zhang 
et al. and this may be due to the fact that we have used more datasets. However, the 
priors of the parameters may also alter the result. A notable difference between the models 
is due to the fact that one of the parameters they have adopted, the optical depth without 
DM decay, is ill-defined in our model. Furthermore in their model, the cutoff of DM decay 
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FIG. 4: The marginalized posterior likelihood of the total optical depth and that of the DM decay 
reionization parameter. 




FIG. 5: The marginalized 2D likelihood contours. The inner and outer contours correspond to 
68.3% and 95.4% confidence levels, respectively. 



reionization at z = 7 was explicitly chosen. These differences might have led to Fig. 1 in 
their paper which shows a preference for a zero dark matter decay rate, a feature that is 
absent in our results. 

The constraint we have obtained is / F < 0- 59 x 10~^^s^^ at 95.4% confidence level. This 
result is about a factor of 3 better than Ref . 1^ . In the limit where reionization is negligible. 



Ref. 18(1 cannot give a strong bound on F because they have ignored the impact of DM decay 
on cosmological perturbations. Hence, their constraint on DM decay is essentially only 
from CMB polarization data. Our analysis, which combines many different cosmological 
datasets and includes the calculation of the impact of DM decay on all the observables, 
gives a stringent constraint on F even in the / = limit. As for the limit of significant 



reionization, our bounds are, as mentioned earlier, a significant improvement over Ref. [18 



The difference may be due to the fact that we have used more datasets. However, the priors 
of the parameters may also alter the result. A notable difference between the models is due 
to the fact that one of the parameters they have adopted, the optical depth without DM 
decay, is ill-defined in our model. Also in their model, the cutoff of DM decay reionization 
at 2; = 7 was explicitly chosen. These differences might have led to Fig. 1 in their paper 
which shows a preference for a zero dark matter decay rate, which is absent in our analysis. 
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4. IMPLICATIONS FOR PARTICLE PHYSICS MODELS WITH DECAYING 
COLD DARK MATTER 



Our results most certainly impose constraints on extensions of the Standard Model of 
particle physics (SM) with DM candidates. Making the assumption that DM decays into 
SM fields, wc will investigate all the probable decay channels unless forbidden by either 
symmetry or kinematics, or highly suppressed by phase space considerations. We then sum 
up their partial decay widths to obtain the functional form for the lifetime of each of the 
decaying DM candidates. One might be worried about including even the decays to non- 
rclativistic particles as that might invalidate our earlier assumption that cosmologically, the 
dark matter decay products are relativistic. While it is true that the DM particle can and 
will decay (provided it is not kinematically forbidden) into non-relativistic massive gauge 
bosons or heavy quarks, these heavy particles themselves will be assumed to subsequently 
decay very rapidly into much lighter particles of the SM that will be relativistic. Obviously 
in specific models, certain channels could be expressly forbidden by symmetries and this can 
also be handled by our analysis. 

To obtain the functional form for the lifetime of the decaying DM, we will approach it 
from the point of view of effective field theory. We will do a model independent analysis 
by considering generic Lagrangian terms for these decays with the corresponding coupling 
constants acting as Wilson coefficients. Since the DM is electrically neutral, the total charge 
of decay products should also be zero. Moreover, the decay rates for different channels are 
dependent on the intrinsic spin of the DM because of possible spin-dependent couplings. 
Below we discuss the decays of DM with spin 0, 1/2 and 1. We will only consider decay 
processes of the lowest order, as higher-order processes involve more vertex insertions and so 
are assumed to be suppressed. Additionally in our approach, we will work in the framework 
where all the gauge symmetries (including those of Grand Unified Theories if present) except 
for SU{3)c X U{l)em are broken and the effects encoded in the couphng constants of the 
effective operators. This can potentially give rise to naturally very small coupling constants 
as they could contain loop factors or powers of very small dimensionless ratios. This is a more 
cost-effective and model-independent way of taking into account the myriad possibilities of 
UV-completing the Standard Model of particle physics. 

Having obtained the functional form of the lifetime in terms of the fundamental parame- 
ters of the underlying particle physics models, we can then compare it it with the numerical 
value obtained from the cosmological analysis of the previous section. This would allow us to 
place definitive bounds on the fundamental parameters of candidate models for the particle 
physics theory beyond the Standard Model. It should be noted that we will be using the 
most conservative 95.4% confidence level bound on the lifetime of the decaying dark matter, 
i.e. without significant reionization. To assume otherwise would require a more complete 
knowledge of the ionization history of the universe than is currently understood. 

Let us now proceed to the case of a generic scalar DM candidate and see how the above 
ideas are implemented. 

4.1. Spin-0 Dark Matter 

We first consider a spin-zero DM candidate, S. Decays into right-handed neutrinos and 
left-handed anti-neutrinos are not considered as the former may be more massive than S. 
Even if they are light enough for S to decay into, we expect the decay into SM (anti-) 
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neutrinos to be dominant. For this decay channel, we can proceed in the same way as in 
the case of neutral pion decay. In the SM, (anti-) neutrinos couple to other matter in the 
form of the chiral current vl'^^vl since S carries no Lorentz index. The lowest dimension 
operator responsible for this decay will be of the form gof'j^^l + ra'^^)fd^S/K, where we 
parameterize as the coupling constant of dimension zero, and A is some cutoff scale. The 
presence or absence of 75 in the operator depends on whether S couples to the SM neutrinos 
in a vectorial or axial- vectorial way. The corresponding decay rate is given by 

where m/ is the mass of the decay product. Here we have assumed that (anti-) neutrinos 
have Dirac mass. 

We focus our attention on DM with mass ms ^ rrif. Then for decay products such 
as (anti-) neutrinos (or other light SM particles), it is safe to make the approximation 
1 — 4mj/m| ~ 1. Since the neutrino is left handed, we take = — 1 . If decays 
dominantly into Vei'ei our lower bound on F then constrains the following parameter, 

— < 1.3 X lO^^^GeV (95.4% confidence level), (20) 

where we have used IGyr^^ = 2.087 x lO^'^^GeV. 

Here we can see that helicity suppression at work. When the mass of the decaying particles 
is very small, the decay of the spin-zero DM candidate will be suppressed as expected. The 
presence of hehcity suppression gives us a value of go that is larger than in most other cases, 
as we will see. For example, if the mass of the DM candidate is m5 ~ lOOGeV, the neutrino 
mass around ~ 2eV, and the cutoff is A ~ lOTeV, then the coupling constant qq has to 
be ~ 10~^^. If on the other hand the coupling constant Qq is 0{1) and we take the same 
values of the neutrino masses and of the DM candidate, then we have that A ~ lO^^GeV. 

Table HI] in the Appendix lists out possible Lagrangian terms for the decay of S into SM 
particles, and the corresponding decay rates, summing over final state spins. Apart from 
focusing solely on the S ^ vv channel, we can consider scenarios in which all the interaction 
terms in Table HTl are present to contribute to the decay rate of S", with all coupling constants 
of the same order of magnitude, Qq. For simplicity, we will also assume rridecay productions is 
negligible compared to 1. This would immediately imply that the helicity-suppressed term, 
5'o/7'*(l + '^a75)/c^Aj5'/A, gives rise to insignificant decay rate when compared to other terms. 
So the most relevant terms are the ones that come from the operator g^S f{l + irp'^^)f . Then 
we have 

/e SM ^ /G SM 

= + (21) 

Here the decay to ti is not included as this channel may not be kinematically feasible for a 
DM particle of ~ lOOGeV. The parameter is also assumed to be the same for all /. In a 
similar way, rates for the other decays into SM gauge bosons can be worked out: 

F(^ ^ 77) + F(5 ^ gg) + V{S ZZ) + F(5 ^ W+W-) + V{S ^ Z7) 
~ (80 + 640 + 82 + 81 + 10). (22) 
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Our result imposes an upper bound on the sum of decay rates via all the channels. Taking 
= as an example, our result would give us the constraint 

glnis < 3.9 x lO^^^GeV (95.4% confidence level). (23) 



4.1.1- Messenger number violation in gauge- mediated supersymmetry breaking theories 

Let us now investigate the messenger parity in gauge-mediated supersymmetry break- 
ing theories. In these models, there could potentially be a dark matter candidate coming 
from the electromagnetically-neutral scalar field that is formed from the SU (2) doublets of 
the 5 and 5 of the messenger sector [68]. However, it is usually not easy to realize this 
because the lightest odd-messenger parity particle (LOMPP) often turns out not to be the 
electromagnetically-neutral field that we require. It has been claimed in the same paper that 
certain F-terms would lift the degeneracy. If we further assume that it does not significantly 
modify the effective low energy theory, the analysis becomes very much model-independent 
as there are only certain couplings that lead to decay of the LOMPP. Following [68|, the 
Kahler potential is given by 

K = Jd^e (51,5m + 51,5m + 51,5^ + lOllOi.) + ^ (sl^lOl + 5l,5,.10^ + /i.e.) , (24) 

where 5m and 5m are the messengers and 5f and lOp are the ordinary superfields. The 
terms that are Planck suppressed are the ones that violate messenger number by one unit. 
As for the superpotential, we have 

W = J d^e p5m5m + jf^^MlO'p, (25) 

where p is the supersymmetry breaking spurion field and once again the terms that are 
Planck suppressed are dimension-5 messenger number violating terms. Without full knowl- 
edge of the UV-sensitive physics (F-terms that lift the other fields while retaining a viable 
LOMPP), we can still give an estimate of the order of magnitude of the decay rate of the 
LOMPP, 

2 3 

r ~ AT^o^^^^Ffc, (26) 

where N are the different degrees of freedom that the LOMPP can decay into. is a 
function that contains the kinematic information and we will assume that is close to one. 
We can then put a constraint on the coupling constant and on the messenger mass. Since 
the lifetime is lOOGyr, then for F^. ~ 1 and ~ 100, we have 

9l {^%r) ^ 6 X 10-^^GeV, (27) 

where we have assumed one universal coupling constant qq. For the sake of discussion, if we 
consider a coupling constant of order one, we get a small messenger mass rrimess ~ 0.02GeV. 
This can be improved if we go to dimension six operators which gives a generic decay rate 
of the form 
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(28) 



Note that instead of a m^g^^/M^ suppression, we now have m^g^^/Mp. This gives a viable 
scenario since the messenger mass now needed is rrimess ~ 4000 TeV, for a couphng constant 
of order one. 



4.2. Spin-1/2 Dark Matter 

We now consider a massive DM of spin-1/2 (let us call it ip) that decays into SM particles. 
Without a specific model, we assume if) decays dominantly via two-body decays and focus 
our attention to this phenomenon. This means must decay into one SM fermion / and 
one SM gauge boson G. Since the DM candidate must be neutral, the posibilities for a 
two body decay of spin 1/2 DM into SM particles arc (/, G) = (z/, Z), (/^, l^"^). For the 
case that ip \s a, Dirac fermion, the two-body decays are mediated by the effective operator 
guGfjifl^i^ + + gh^^ii^l^ + ^75)/- The first term gives rise to ip decay while the 

second one is responsible for the decay of -0. Again summing over the final state spins and 
averaging over the spin of the decaying we find the decay rate to be 

r(^ ^ }G) = r(^ - JG) = ^I^Jxl"-^,^) L (^,^)+ r'. ( ^, 

(29) 

where A(a,6) = {I + a -h){l - a -h){l - a + h){l + a + h) anda;(a,6) = (l + a-6)(l-a- 
6)[2a^ + (1 + bY]) and denotes respectively the mass of particle A. 

Now consider the case where the decay of the fermionic DM candidate comes from an 
operator guipH{l + irpj5)f + g})fH{l + irp'y5)ip, where H is the Higgs boson^ of the SM. 
In this case the decay rate of ^0 — > if -|- 1/ is given by 



r(^ ^ H.) = r(^ ^ HV) = I V rn,\ 

167r y m^J \_ J ^ \ J 

(30) 

where z{a, h) ^ I -\- - h'^ -\- 2a. 

We can now consider a simple scenario in which r = rp = for all the decay channels 
and they all have the same coupling constant go- Then the total decay rate of iJj is given by 
summing over all the possible channels: 

r(V' ^ 2 body) = 3 r(V' ^ Zu) + r(V' ^ w+e-) + r(V' ^ w+i^-) 

+r(v^ ^ W+T-) + 3 r(V' ^ Hu) (31) 

= 126.5 i^^r 



We consider the Higgs in this particular case because this interaction would arise at the same or lower 
order than the other one we considered. With decaying DM of other spins, there would be an additional 
suppression from the ratio of electroweak scale over the cut-off scale. 
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where we have picked ~ 200GeV, H ~ lOOGeV and assumed all three generations of 
neutrinos have masses ~ leV. The factor 3 in Equation (13T|) is for three generations of 
neutrinos. Our cosmological bound then gives us the constraint 

Dirac fermion : {gol < 4.0 x 10"^^ (95.4% confidence level). (32) 

For the case that ip is a Majorana fermion, the above analysis follows through. Given the 
same interaction terms as those shown above, the partial decay rates of a Majorana ip are 
exactly the same as Equations (!29|) and (!30l) . There are, however, no distinction between ijj 
and ip in this case any more. In a four-component spinor notation, ip and ip relate to each 
other via the charge conjugation matrix. This means the total decay rate of a Majorana ip 
has contributions from decays into 'particles' and decays into 'anti-particles'. 

In the same simple scenario we considered above, the total decay rate of a Majorana ip 
will be increased by a factor of 2 compared to Equation fl3T|) . The constraint on the coupling 
constant will correspondingly be tightened by a factor of 

Majorana fermion : IqmI < 2.8 X 10"^^ (95.4% confidence level). (33) 



4-2.1. R-parity violation in minimal supergravity models 

Undoubtedly, the most thoroughly investigated models in the supersymmetric menagerie 



are the minimal supergravity (mSUGRA) models [69|, l70|, lDJ, l72|, |73| . While the theoretical 
motivation for universality of scalar masses, gaugino masses and trilinear terms is ques- 
tionable (since these values depend on the mechanism by which supersymmetry breaking is 
transmitted to our sector), it has nevertheless remained a useful benchmark. For our pur- 
poses, it is sufficient for us to use the fact that in a variety of these mSUGRA models, the 
lightest supersymmetric particle (LSP) is a neutral particle that is overwhelmingly composed 
of the spin-1/2 supersymmetric partner of the B- gaug e boson called the bino, B. There are 
of course technically natural classes of models [7J, |75|, very similar to mSUGRA theories 
that will give bino as the LSP, and the analysis below would similarly apply to them. 

In the presence of R-parity violation, the bino LSP would of course decay. Traditionally, 
theories with R-parity violation were often assumed to be unable to provide a dark matter 
candidate. Here, we can turn this around and ask what the couplings of the theory have to 
be so that the theory can still furnish us with viable dark matter candidate. To do that, we 
need to first explore the possible decays. 

While the two-body decay might seem to have a more favorable phase space, these decays 
however would arise from Feynman diagrams [t^I only if we have the R-parity violating terms 
together with the introduction of an additional loop and further suppression by dimensionless 
ratios of electroweak scale over the cut-off scale. We will therefore assume that the bino 
will dominantly decay into three SM particles via the trilinear R-parity violating terms. If 
for some particular models, one needs to add in some of the two-body decay terms, one can 
look up the Appendix or the previous subsection for the relevant cross-sections and include 
them in the overall analysis. 

Neglecting all final state masses, the decay rate for a three-body bino decay is given by 

/ dE^ / dE2 Y,\M\^ (34) 

H^B-^^ spins 
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where Ei is the energy of the final particle i^, and the summation symbol means averaging 
over initial spins and summing over final spins. The amplitudes squared for three-body 
decays of neutralino due to trilinear R-parity violating terms have been evaluated and shown 



m 



78|, |79(], with the appropriate spin summing/averaging. Strictly speaking a neutralino 
is a superposition of the bino and three other fermionic supersymmetric particles but for 
our purposes, it is sufficient to consider the bino to the lightest neutralino and the LSP. The 
results in ^79| can be easily applied to LSP decay by demanding the neutralino has a 100% 



bino component, i.e. by setting A^^i = 1 and N^n = 0, = 2, 3, 4 in the notation of ^]. For 
simplicity, all final state masses are neglected in the analysis below. We have also ignored 
the mixings and the widths of the sfermions, which mediate the decay as internal lines in 
the Feynman diagrams. 

Given the R-parity violating superpotential term 



W, 



LLE 



(35) 



(where i,j and k, each of which runs from 1 to 3, are generation indices, a and p_are SU{2)l 
indices, and the superscript c indicates charge conjugation), the decay channel B ^t^j^k 
is possible. Using the generic expression for amplitude squared in and putting in our 
simplifications, we get the decay rate 
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(36) 



where g' is the gauge coupling of U{1)y, mj^ is the mass of the scalar superpartner of particle 

fn and Ys denotes the hypercharge of a superfield S (for example. Ye = —1). K{x) and 
P{x,y) are functions defined in the Appendix. 

Another trilinear R-parity violating superpotential term is 



(37) 



with the SU{3)c index a. This term gives rise to the decays B efujdk and B —>■ 
Vidjdk- The decay rates for these channels are similar to the one above, with the appropriate 
substitution of superpartner masses and prefactors: 
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(38) 



^ Of course the identification of a final particle as particle i is arbitrary. This arbitrariness does not change 
the final expression for F. 
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and 
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(39) 



Note that the numerical value of an SU (3)c colour factor has been included in the prefactors 
of Equations ([38D and ([39D. 

In a similar way, the decay channel B Uidjdk is allowed by the superpotential term 



where a, jd and 7 are all SU{?>)c indices. The corresponding decay rate is 

48 



(40) 
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(41) 



Here a different SU{'i)c colour factor has been included in the prefactor. 

It should be pointed out that the bino is a Majorana fermion. This means what we shown 
above is only half of its possible decay channels: the other decay channels are obtained by 
applying charge conjugation to all the final particles in any of the above channels. The 
decay rates, however, are invariant under charge conjugation. 

Because the LLE term contains two copies of L's and they contract with the Levi-Civita 
tensor, Ajj^ is anti-symmetric in i and j. Thus it only represents nine couplings. Similarly, 
X'l^k is anti-symmetric in j and k. This argument is not applicable to X'^.^i so it does indeed 



rn^ — I I — h 

contain 27 couplings (see, for example, [6l. I78l. |79[|). 

As a simple application of our cosmological constraint on the DM decay rate, we assume 
~ lOOGeV and all the sfermions have masses ~ SOOGeV. We also assume all the 
non-zero R-parity violating couplings attain the same value A, i.e. 



A,: 



i2j2k2 



(42) 



Summing over all the possible 3-body decay channels of bino, the total decay rate is given 
by 



T{B 3 body) = 2[9T{B ^ e+z/^e^ ) + 27 T{B eju^dk) 
+27 r{B Uidjdk) + 9 r{B uidjdk) ] 
= 0.00144|A|2 



(43) 



where we have used g' = 0.36. Our cosmological bound then constrain the coupling constant 
to be 

|A| < 1.2 X 10"^° (95.4% confidence level). (44) 
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In comparison, one of the strongest constraint on R-parity violation comes from the 
consequent baryon number violation that arises due to the former. Ref . [80| gives a value of 
A" < 10~^ for the most constrained of all the A's. So if indeed the assumption that we have 
bino-like DM holds true, then the most stringent limits on R-parity violation would come 
from our analysis. 



4.3. Spin-1 Dark Matter 

We now consider a massive DM of spin-1 (let us call it x) that decays into SM particles. 
Since the x'^ fi^ld carries one Lorentz index, it contracts with other SM fields differently 
from the spin-0 DM, thus giving rise to different interaction terms and decay rates. 

In contrast to the decay of spin-0 DM, helicity suppression is not observed in the decay 
of X — > uu. x^ can directly coupled to the neutrino current I'Llfi^L, without any insertion of 
d^. On the other hand, every spin-1 particle has to obey the Landau- Yang theorem [sil, 8^ 



which states that because of rotational invariance, it cannot decay into two massless spin-1 
particles. Hence, the decays x ^ 11 and x ^ 99 a^re not allowed. The possible partial 
decay widths (with summing over final state spins and averaging over the initial state spin) 
for a spin-1 DM are rather numerous and not that illuminating to list them all here. So 
we have relegated them to Table |V] in the Appendix. In the case where the DM is indeed 
an additional U(l) gauge field that is massive, the possibility of kinetic mixing with the 
photon [ssj must be considered. Such a term could be radiatively generated via exchange 



of a field that is charged under both U(l)'s. Following Refs. |84j. l85l| . we can manipulate the 
Lagrangian into a form where the mixing manifests itself in the coefficients of the following 
terms, 9iXiifl^i)- + ^la)!- But this is a term that has already been considered in Table IVl 
To get a feel for the numbers involved, let us now consider a simple model where all the 
interaction terms in Table IVl exist, with all the coupling constants real and of the same order 
of magnitude. Again, we will also assume rudecay product/ Then 

T{x Z^) + T{x - ZZ) + T{x - W^W-) + \ + (45) 

9o7r V mi mi mf 



Because rrif/m^ is small for the value of we are considering. 
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An 



Similar to the case of spin-0 DM, the decay to tt is not included here, and the parameter r 
is also assumed to be the same for all /. Note also that we have different m^ dependence 
for the decays into fermion-antifermion and massive gauge bosons, unlike in the spin-0 case. 

Our result then gives an upper bound on the sum of all the decay rates into SM particles. 
For illustration, we consider m^ ~ lOOGeV and r = 0. Our bound on F can then be 
translated into a constraint on the coupling constant: 

gi < 5.8 X 10"^^ (95.4% confidence level), (47) 

where we have substituted mw = 80GeV and mz = 91GeV. 
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4-3.1. T-Parity violation in little Higgs models 



Little Higgs models with T-parity violation is another possible scenario in which the dark 
matter candidate decays. Analogous to R-parity in SUSY models, all non-SM particles in 
Little Higgs model are assigned to be T-odd, while all SM ones T-even. The T-parity then 
requires all coupling terms to have an even number of non-SM fields. This forbids the contri- 
bution of the non-SM particles to the oblique electroweak parameters, and consequently, the 
symmetry breaking scale / can be lowered to about ITeV |86!]. The Lightest T-odd Particle 
(LTOP), moreover, is stable and has often been nominated as a dark matter candidate. 
However, Ref. [sj has pointed out that anomalies in general give rise to a Wess-Zumino- 

have constructed Little Higgs 



Witten (WZW) term, which breaks the T-parity (Refs.[87. 
models free of the usual WZW term). This means the LTOP is not exactly stable. Indeed, 
phenomenological consequences of the WZW term in the Littlest Higgs model have been 
studied in Ref.fSQ. In their model, the LTOP is the massive partner of photons (denoted 
by Ah) and the WZW term contains direct couplings of Ah to the Standard Model W 



bosons, Z bosons and photons. Ref.[90|], moreover, pointed out that such couplings can 
generate two-body decay of Ah to SM fermions. Ah — //, via triangular loop diagrams. 

In an attempt to be as model-independent as possible, we parameterize the couplings of 
Ah to the SM gauge bosons as 



(48) 

where / is the symmetry breaking scale, g' the U{1) gauge coupling, Nz, Ny/ and N^z are 
numbers whose values depend on the exact realization and the UV completion. 

Generically, the mass of Ah is proportional to /. If we take / to be the natural symmetry 
breaking scale (i.e. ~ ITeV) in Little Higgs models, then m^^ > 2mz. As an example, in 
(sol . 90], we have 

246GeV is the Higgs vev. The condition for > 2mz is satisfied when 
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(49) 



where v 
g' ~ 0.36 and / > 1165GeV. 

In the case of > 2m z, the decay channels of Ah ZZ and Ah — > W^W~ 

are kinematically allowed, and, for simplicity, we assume these processes (together with 
Ah — >■ Z'^) to be the dominant ones. With the interaction terms in equation (HHj) . the decay 
rates for these channels at the lowest order are given by 
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(50) 
(51) 
(52) 



The sum of these decay rates is then constrained by our bound on the dark matter lifetime. 
For the Littlest Higgs model, N^z = and Nw = Nz- The sum of the above decay rates is 
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then reduced to 

r = (53) 

where we have used the approximations m^y ~ m^, vhah — g'f/\/5 and have neglected all 
the mass ratios. Our bound on F then gives us the constraint 

< 4.7 X 10-^^GeV (95.4% confidence level), (54) 

which of course is not reasonable as we typically expect Nz ~ 1 and / ~ 1 TeV. But it 
vividly illustrates the utility of our approach when it comes to ruling out particle physics 
models that claim to have dark matter candidates. 



4.4. General Dimensional Considerations 

We can draw some generalizations from the above cases if wc do a simple dimensional 
analysis. The coupling of a spin-0, spin-1/2 or spin-1 dark matter candidate S to an operator 
O can be parameterized, with suppression of indices and 0{1) factors, as 

L =) ^X^O, (55) 

where g is a dimensionless coupling constant, A is the scale where unknown new physics is 
integrated out to give the operator O, and n is the sum of the dimensions of S and O. The 
decay rate for such dark matter candidate is given in general by 

r - J^s^f''F„ (56) 

where Fk is a function that contains the kinematics of the decay, assumed to be of order 
one for simplicity. For n — 5, ms — lOOGeV and g ~ C(l), the cutoff scale should be of the 
order of A > lO^^GeV, suggesting that we must go to operators of higher dimensions and 
thus more A suppression. For n = 6, ms = lOOGcV and g ~ the cutoff scale can be 

as low as A ~ lO^^GeV. On the other hand, for n = 5, if the cutoff is taken at the Planck 
scale (A ~ lO^^GeV) and we keep the same value of m^, the couphng constant can only be 
as large as 51 ~ 10~^. Finally to recover the cases discussed above for spin-0, spin-1/2 and 
spin-1 particles, we can take n — A and ms ~ 100 GeV to give us a coupling constant as 
large as g ^ 10~^^. 

A few words should be reiterated about the smallness of the coupling constant. We had 
taken an extremely conservative value for our cutoffs, usually ~ 10 TeV. In an effective 
theory with a low cutoff arising from a high scale fundamental theory, say at Planck scale, 
there will be a multitude of effective operators containing mass insertions (leading to small 
dimensionless ratios such as m/Mp) or loops (giving factors of I/IGtt^), making these tiny 
coupling constants natural. The small dimensionless ratios could arise from, say, the decay 
being mediated by some massive field much like what we have in proton decay via exchange 
of heavy X bosons in the context of Grand Unified Theories. The onus is then on the model 
builders to refine their models in a technically natural way to satisfy the constraints we have 
derived above without having to compromise other phenomenological constraints on their 
models. 
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5. CONCLUSIONS AND OUTLOOK 



We have performed a full cosmological analysis using the available datasets from cosmic 
microwave background, Type la supernova, Lyman-a forest, galaxy clustering and weak 
lensing observations to determine the extent by which we can constrain decaying dark matter 
models which are very typical in most extensions of the Standard Model of particle physics. 

In the scenario where there is negligible reionization of the baryonic gas by the decay- 
ing dark matter, we have found that the late-time Integrated Sachs- Wolfe effect gives the 
strongest constraint. The lifetime of a decaying dark matter has the bound F"^ > lOOGyr 
(at 95.4% confidence level). Because of cosmic variance, the results are not likely to improve 
significantly with the WMAP-9yr data. 

When there is significant reionization of the baryonic gas due to the decaying dark matter, 
the bounds become more restrictive as the CMB polarization is well measured. In this 
scenario, the lifetime of a decaying dark matter is (/F)"^ > 5.3 x 10^ Gyr (at 95.4% 
confidence level) where / is a phenomenological factor related to the degree of reionization. 
With even more CMB polarization data, one could conceivably distinguish the reionization 
due to decaying dark matter from reionization due to star formation, thereby giving us even 
better bounds on the lifetime of the dark matter. We expect that the the 21cm cosmological 
observation in the future would give us even greater precision as it is expected to probe the 
reionization history at redshifts 6 < z < 30. 

Having obtained the cosmological constraints, we turned our attentions to the particle 
physics aspects of it. For completeness and motivated by the utility of such an exercise, 
we systematically tabulated the decay cross- sections for a spin-0, spin-1/2 and spin-1 dark 
matter candidate into the Standard Model degrees of freedom. This enabled us to simply 
sum up all the relevant contributions for a particular model of particle physics and arrive at 
the functional form of the hfetime of the decaying dark matter. We repeated this process 
for a variety of representative models from the following classes of theories: generic super- 
sysmmctric scenario, gauge-mediated supersymmetry breaking models and the little Higgs 
theories. Imposing the limits from our cosmological analysis, wc find that gcncrically for 
most models we have looked at, the dimensionless coupling for a decaying dark matter to 
Standard Model fields should be smaller than 10~^^. 

This restriction can be slightly relaxed if the dark matter decays solely into light particles 
via helicity suppressed interaction terms, in which case, the small mass of the decay prod- 
ucts suppresses the decay rate. If, for instance, the dark matter decays purely via hehcity 
suppressed terms into uP with Dirac mass of ~ 2eV, then the dimensionless coupling can be 
as large as 10~^^. In addition to constraining the coupling, one can assume it to be of 0{1) 
and estimate the scale of new physics which suppresses the decay rate. In all cases, either 
the coupling attains a small value or the new physics come from a huge scale, both of which 
would need interesting and exotic physics to realize if indeed the dark matter does decay 
via dimension-4 or dimension-5 terms. In the case of exclusive helicity suppressed decays, 
moreover, one has to explain why other interaction terms are absent in the model. A more 
promising avenue, which we briefly mentioned in the previous section, is to look at models 
where the dark matter decays via dimension-6 operators. The Large Hadron Collider might 
provide us with the identity for dark matter in the very near future, but on the basis of our 
analysis, there will still be much to understand about physics of the dark matter sector and 
how it interacts with the Standard Model. 

In a future paper, we hope to address some of the astrophysical issues of decaying dark 
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matter. The recent spate of results from astrophysical experiments 91, ol, 93, 94] has given 
us much to ponder. The immediate goal would of course be to combine all the astrophysical 
datasets with the cosmological ones that we have considered in the present paper and arrive 
at a set of characteristics that a phenomenologically viable decaying dark matter must 
possess. However, we expect considerable tension between the two classes of constraints. 
The astrophysical ones require that the decays to be significant enough to account for the 
as-yet-unexplained phenomena, while the cosmological ones need decays to be small enough 
because of the late time ISW effect and the CMB polarization observations. To reconcile 
and resolve these two seemingly conflicting classes of observations could be the defining 
challenge of dark matter physics in the next decade. 
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APPENDIX: COMPENDIUM OF DECAY RATES 

This appendix summarizes the lowest order decay rates due to various generic interaction 
terms in the Lagrangian, averaging over the spin of the DM and summing over the spins 
of the decay products. Tables [III HIl and |V] respectively tabulate the decay of a spin-0 DM 
particle S, a spin-1/2 DM particle and a spin-1 DM particle into SM particles. In a 
model-independent way, we write down generic Lagrangian terms which describe possible 
decay channels of DM particles to SM particles. The exact mechanisms which mediate these 
decays are captured by the dimensionless coupling constants, go, go and gi. The reality 
of the interaction terms requires the coupling constants to be real, except in the case of 
X W^W~ and the decay of ip, in which a complex coupling constant is possible. 

We follow standard conventions in denoting our fields. represents a gauge field while 
Y^^ is the corresponding field strength tensor. For various decay channels, Nf and Ng 
respectively denote the number of colours of a fermion species / and a gluon g. r, rp and 
are parameters that describe the relative size of two interaction terms. 



TABLE II: Decay Rate of Spin-0 DM via Different Interaction Terms 
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TABLE III: Two-Body Decay Rate of Spin-1/2 DM via Generic Interaction Terms 



Interaction Term 


Decay Rate 




r(V, ^ fG) = Tii; ^ 

_ \9D\'^ml 


fG) 


+9DfH{l + irp-i^)ip 


r(v^ ^ Hf) = Tii^ ^ 

167r Y 
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TABLE IV: Functions used for the Analysis of Bino Decay 



K{x) 
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2 I 2 



6 x2 



+ 2^ In(x2 + y2_i) [_1 ln(a;2+y2_i)+ ln(x2-l)] 



TABLE V: Decay Rate of Spin-1 DM via Different Interaction Terms 



Interaction Term 


Decay Rate 


5iXm/7^(1 + ^75)/ 


r(x 




- 127r ^X\l ^ X 


"l + 24 + r2 fl-44Y 






r(Y 


- //) 


'A 








77 or gg) 


= 

forbidden by the Landau- Yang theorem 




r(x 


ZZ) 






gief.^p^X^Z-d^ZP 


r(x 


ZZ) 


ybTrm^ \^ mj y 




giW+W-'^d.x^ 

+giw-w+'^d,x^ 


r(x^ 


W+W-) 


1927rm^ V J 

X {4[i?e(<7i)F^ + [/m(rf (l + 4^)} 




r(x^ 


W+W-) 


4^ 

487rm^ Y "^x 




+g*i^,..pax^w-^d''w+p 






X { mg,)]^ (1 - 4^)V [/m(,0]^ (1 + 2^) } 




r(\ 




= Ami. (1 + 41 f 1 




gie^upaX'^ZT'^P 


r(x 


^^7) 







[1] G. Hinshaw, J. L. Wciland, R. S. Hill, N. Odcgard, D. Larson, C. L. Bennett, J. Dunkley, 

B. Gold, M. R. Grcason, N. Jarosik, et al., 803 (2008), 0803.0732. 
[2] E. Komatsu, J. Dunkley, M. R. Nolta, C. L. Bennett, B. Gold, G. Hinshaw, N. Jarosik, 

D. Larson, M. Limon, L. Page, et al., 803 (2008), 0803.0547. 
[3] For a review, see M. Schmaltz and D. Tucker-Smith, Ann. Rev. Nucl. Part. Sci. 55, 229 (2005), 

hep-ph/0502182. 
[4] For a review, see S. P. Martin (1997), hep-ph/9709356. 
[5] C. T. Hill and R. J. Hill, Phys. Rev. D76, 115014 (2007), 0705.0697. 
[6] For a review, see R. Barbier et al., Phys. Rept. 420, 1 (2005), hep-ph/0406039. 
[7] W. Israel, Phys. Rev. 164, 1776 (1967). 
[8] W. Israel, Commun. Math. Phys. 8, 245 (1968). 
[9] B. Carter, Phys. Rev. Lett. 26, 331 (1971). 
[10] R. M. Wald, Phys. Rev. Lett. 26, 1653 (1971). 
[11] T. Banks, Nucl. Phys. B323, 90 (1989). 

[12] L. M. Krauss and F. Wilczek, Phys. Rev. Lett. 62, 1221 (1989). 

[13] R. Cen (2000), astro-ph/0005206. 

[14] F. Borzumati, T. Bringmann, and P. Ullio, Phys. Rev. D77, 063514 (2008), hcp-ph/0701007. 
[15] K. Ichiki, M. Oguri, and K. Takahashi, Phys. Rev. Lett. 93, 071302 (2004), astro-ph/0403164. 
[16] Y. Gong and X. Chen (2008), 0802.2296. 
[17] M. Lattanzi, AIP Conf. Proc. 966, 163 (2008). 

[18] L. Zhang, X. Chen, M. Kamionkowski, Z.-g. Si, and Z. Zheng, Phys. Rev. D76, 061301 (2007), 
0704.2444. 

[19] S. Palomares-Ruiz, Phys. Lett. B665, 50 (2008), 0712.1937. 

[20] H. Yuksel and M. D. Kistler, Phys. Rev. D78, 023502 (2008), 0711.2906. 

[21] M. J. Strassler and K. M. Zurek, Phys. Lett. B651, 374 (2007), hep-ph/0604261. 

[22] C.-P. Ma and E. Bertschinger, Astrophys. J. 455, 7 (1995), astro-ph/9506072. 

[23] L. Kofman, D. Pogosyan, and A. A. Starobinsky, Sov. Astron. Lett. 12, 175 (1986). 

[24] J. R. Bond and G. Efstathiou, ApJ. Lett. 285, L45 (1984). 

[25] X.-L. Chen and M. Kamionkowski, Phys. Rev. D70, 043502 (2004), astro-ph/03 10473. 

[26] D. R. G. Schleicher, R. Banerjee, and R. S. Klesser, Phys. Rev. D78, 083005 (2008), 0807.3802. 

[27] S. Seager, D. D. Sasselov, and D. Scott, ApJ. Lett. 523, LI (1999), arXiv:astro-ph/9909275. 

[28] S. Seager, D. D. Sasselov, and D. Scott, ApJ. Suppl. 128, 407 (2000), arXiv:astro-ph/9912182. 

[29] A. Lewis and S. Bridle, Phys. Rev. D 66, 103511 (2002), arXiv:astro-ph/0205436. 

[30] W. C. Jones, P. A. R. Ade, J. J. Bock, J. R. Bond, J. Borrih, A. Boscaleri, P. Cabella, C. R. 

Contaldi, B. P. Crill, P. de Bernardis, et al.. New Astronomy Review 50, 945 (2006). 
[31] F. Piacentini, P. A. R. Ade, J. J. Bock, J. R. Bond, J. Borrill, A. Boscaleri, P. Cabella, C. R. 

Contaldi, B. P. Crih, P. de Bernardis, et al., Astrophys. J. 647, 833 (2006), arXiv:astro- 

ph/0507507. 

[32] T. E. Montroy, P. A. R. Ade, J. J. Bock, J. R. Bond, J. Borrill, A. Boscaleri, P. Cabella, C. R. 

Contaldi, B. P. Crih, P. de Bernardis, et al., Astrophys. J. 647, 813 (2006), arXiv:astro- 

ph/0507514. 
[33] C. L. Reichardt et al. (2008), 0801.1491. 

[34] C. L. Kuo, P. A. R. Ade, J. J. Bock, J. R. Bond, C. R. Contaldi, M. D. Daub, J. H. Goldstein, 



26 



W. L. Holzapfel, A. E. Lange, M. Lueker, et al. (2006), astro-ph/0611198. 
[35] M. C. Runyan, P. A. R. Adc, R. S. Bhatia, J. J. Bock, M. D. Daub, J. H. Goldstein, C. V. 
Haynes, W. L. Holzapfel, C. L. Kuo, A. E. Lange, et al., ApJ. Suppl. 149, 265 (2003), 
arXiv:astro-pli/0303515. 

[36] J. H. Goldstein, P. A. R. Ade, J. J. Bock, J. R. Bond, G. Cantalupo, G. R. Gontaldi, M. D. 
Daub, W. L. Holzapfel, G. Kuo, A. E. Lange, et al., Astrophys. J. 599, 773 (2003), arXiv:astro- 

ph/0212517. 

[37] T. J. Pearson, B. S. Mason, A. G. S. Readhead, M. G. Shepherd, J. L. Sievers, P. S. Udom- 
prasert, J. K. Gartwright, A. J. Farmer, S. Padin, S. T. Myers, et al., Astrophys. J. 591, 556 

(2003) , arXiv:astro-ph/0205388. 

[38] A. G. S. Readhead, B. S. Mason, G. R. Gontaldi, T. J. Pearson, J. R. Bond, S. T. Myers, 
S. Padin, J. L. Sievers, J. K. Gartwright, M. G. Shepherd, et al., Astrophys. J. 609, 498 

(2004) , arXiv:astro-ph/0402359. 

[39] A. G. S. Readhead, S. T. Myers, T. J. Pearson, J. L. Sievers, B. S. Mason, G. R. Gontaldi, J. R. 
Bond, R. Bustos, P. Altamirano, G. Achermann, et al., Science 306, 836 (2004), arXiv:astro- 
ph/0409569. 

[40] J. L. Sievers, G. Achermann, J. R. Bond, L. Bronfman, R. Bustos, G. R. Gontaldi, G. Dickin- 
son, P. G. Ferreira, M. E. Jones, A. M. Lewis, et al., Astrophys. J. 660, 976 (2007). 

[41] G. Dickinson, R. A. Battye, P. Garreira, K. Gleary, R. D. Davies, R. J. Davis, R. Genova- 
Santos, K. Grainge, G. M. Gutierrez, Y. A. Hafez, et al., MNRAS 353, 732 (2004), 
arXiv:astro-ph/0402498. 

[42] N. W. Halverson, E. M. Leitch, G. Prykc, J. Kovac, J. E. Garlstrom, W. L. Holzapfel, 
M. Dragovan, J. K. Gartwright, B. S. Mason, S. Padin, et al., Astrophys. J. 568, 38 (2002), 
arXiv:astro-ph/0104489. 

[43] E. M. Leitch, J. M. Kovac, N. W. Halverson, J. E. Garlstrom, G. Pryke, and M. W. E. Smith, 
Astrophys. J. 624, 10 (2005), arXiv:astro-ph/0409357. 

[44] S. Hanany, P. Ade, A. Balbi, J. Bock, J. Borrill, A. Boscaleri, P. de Bernardis, P. G. Ferreira, 
V. V. Hristov, A. H. Jaffc, ct al., ApJ. Lett. 545, L5 (2000), arXiv:astro-ph/0005123. 

[45] J. R. Bond, G. R. Gontaldi, U.-L. Pen, D. Pogosyan, S. Prunet, M. I. Ruetalo, J. W. Wadsley, 
P. Zhang, B. S. Mason, S. T. Myers, et al., Astrophys. J. 626, 12 (2005), arXiv:astro- 
ph/0205386. 

[46] M. Kowalski et al. (2008), 0804.4142. 

[47] S. Golc, W. J. Percival, J. A. Peacock, P. Norberg, G. M. Baugh, G. S. Frenk, I. Baldry, 
J. Bland-Hawthorn, T. Bridges, R. Gannon, et al., MNRAS 362, 505 (2005), arXiv:astro- 
ph/0501174. 

[48] M. Tegmark, D. J. Eisenstein, M. A. Strauss, D. H. Weinberg, M. R. Blanton, J. A. Frieman, 
M. Fukugita, J. E. Gunn, A. J. S. Hamilton, G. R. Knapp, et al., Phys. Rev. D 74, 123507 
(2006), arXiv:astro-ph/0608632. 

[49] D. J. Eisenstein, I. Zehavi, D. W. Hogg, R. Scoccimarro, M. R. Blanton, R. G. Nichol, R. Scran- 
ton, H.-J. Seo, M. Tegmark, Z. Zheng, et al., Astrophys. J. 633, 560 (2005), arXiv:astro- 
ph/0501171. 

[50] W. J. Percival, R. G. Nichol, D. J. Eisenstein, J. A. Frieman, M. Fukugita, J. Loveday, A. G. 
Pope, D. P. Schneider, A. S. Szalay, M. Tegmark, et al., Astrophys. J. 657, 645 (2007), 
arXiv:astro-ph/0608636. 

[51] R. Massey, J. Rhodes, A. Leauthaud, P. Gapak, R. Ellis, A. Koekemoer, A. Refregier, N. Scov- 
ille, J. E. Taylor, J. Albert, et al., ApJ. Suppl. 172, 239 (2007), arXiv:astro-ph/0701480. 



27 



[52] J. Lesgourgucs, M. Vicl, M. G. Hachnclt, and R. Massey, Journal of Cosmology and Astro- 
Particle Physics 11, 8 (2007), arXiv:0705.0533. 

[53] H. Hoekstra, Y. Mellier, L. van Waerbeke, E. Semboloni, L. Fu, M. J. Hudson, L. C. Parker, 
I. Tereno, and K. Benabed, Astrophys. J. 647, 116 (2006), arXiv:astro-ph/0511089. 

[54] J. Benjamin et al., arXiv: astro-ph/0703570 (2007), astro-ph/0703570. 

[55] C. Schimd, I. Tereno, J.-P. Uzan, Y. Mellier, L. van Waerbeke, E. Semboloni, H. Hoekstra, 

L. Fu, and A. Riazuelo, A&A 463, 405 (2007), arXiv:astro-ph/0603158. 
[56] M. Hetterscheidt, P. Simon, T. Erben, P. Schneider, M. Schirmer, J. P. Dietrich, H. Hilde- 

brandt, O. Cordes, T. Schrabback, L. Haberzettl, et al.. The Messenger 126, 19 (2006). 
[57] M. Hetterscheidt, P. Simon, M. Schirmer, H. Hildebrandt, T. Schrabback, T. Erben, and 

P. Schneider, A&A 468, 859 (2007), arXiv:astro-ph/0606571. 
[58] H. Hoekstra, H. K. C. Yee, M. D. Gladdcrs, L. F. Barrientos, P. B. Hall, and L. Infante, 

Astrophys. J. 572, 55 (2002), arXiv:astro-ph/0202285. 
[59] H. Hoekstra, H. K. C. Yee, and M. D. Gladders, Astrophys. J. 577, 595 (2002), arXiv:astro- 

ph/0204295. 

[60] L. Van Waerbeke, Y. Mellier, and H. Hoekstra, A&A 429, 75 (2005), arXiv:astro-ph/0406468. 

[61] N. Kaiser, Astrophys. J. 388, 272 (1992). 

[62] N. Kaiser, Astrophys. J. 498, 26 (1998), arXiv:astro-ph/9610120. 

[63] M. Viel, M. G. Haehnelt, and V. Springel, MNRAS 354, 684 (2004), arXiv:astro-ph/0404600. 
[64] T.-S. Kim, M. Viel, M. G. Haehnelt, R. F. Carswell, and S. Cristiani, MNRAS 347, 355 

(2004), arXiv:astro-ph/0308103. 
[65] R. A. C. Croft, D. H. Weinberg, M. Bolte, S. Buries, L. Hernquist, N. Katz, D. Kirkman, and 

D. Tytler, Astrophys. J. 581, 20 (2002), arXiv:astro-ph/0012324. 
[66] P. McDonald, U. Seljak, R. Cen, D. Shih, D. H. Weinberg, S. Buries, D. P. Schneider, D. J. 

Schlegel, N. A. Bahcall, J. W. Briggs, et al., Astrophys. J. 635, 761 (2005), arXiv:astro- 

ph/0407377. 

[67] P. McDonald, U. Seljak, S. Buries, D. J. Schlegel, D. H. Weinberg, R. Cen, D. Shih, J. Schaye, 
D. P. Schneider, N. A. Bahcall, ct al., ApJ. Suppl. 163, 80 (2006), arXiv:astro-ph/0405013. 
[68] S. Dimopoulos, G. F. Giudice, and A. Pomarol, Phys. Lett. B389, 37 (1996), hep-ph/9607225. 
[69] A. H. Chamseddine, R. L. Arnowitt, and P. Nath, Phys. Rev. Lett. 49, 970 (1982). 
[70] R. Barbieri, S. Ferrara, and C. A. Savoy, Phys. Lett. B119, 343 (1982). 
[71] L. J. Hall, J. D. Lykken, and S. Weinberg, Phys. Rev. D27, 2359 (1983). 
[72] E. Crcmmcr, P. Fayet, and L. Girardello, Phys. Lett. B122, 41 (1983). 
[73] N. Ohta, Prog. Theor. Phys. 70, 542 (1983). 

[74] H.-S. Goh, M. A. Luty, and S.-P. Ng, JHEP 01, 040 (2005), hep-th/0309103. 
[75] H.-S. Goh, S.-P. Ng, and N. Okada, JHEP 01, 147 (2006), hep-ph/0511301. 
[76] S.-P. Ng and N. Okada, JHEP 09, 040 (2007), 0705.2258. 
[77] S. Dawson, Nucl. Phys. B261, 297 (1985). 

[78] H. K. Dreiner and P. Morawitz, Nucl. Phys. B428, 31 (1994), hep-ph/9405253. 
[79] E. A. Baltz and P. Gondolo, Phys. Rev. D57, 2969 (1998), hep-ph/9709445. 
[80] M. Chemtob, Prog. Part. Nucl. Phys. 54, 71 (2005), hep-ph/0406029. 
[81] L. D. Landau, Dokl. Akad. Nauk Ser. Fiz. 60, 207 (1948). 
[82] C.-N. Yang, Phys. Rev. 77, 242 (1950). 
[83] B. Holdom, Phys. Lett. B166, 196 (1986). 

[84] K. S. Babu, C. F. Kolda, and J. March-Russell, Phys. Rev. D54, 4635 (1996), hcp-ph/9603212. 
[85] K. S. Babu, C. F. Kolda, and J. March-Russell, Phys. Rev. D57, 6788 (1998), hep-ph/9710441. 



28 



H.-C. Cheng and I. Low, JHEP 08, 061 (2004), hcp-ph/0405243. 

C. Csaki, J. Heinonen, M. Perelstcin, and C. Spcthmann (2008), 0804.0622. 

D. Krohn and I. Yavin, JHEP 06, 092 (2008), 0803.4202. 

V. Barger, W.-Y. Keung, and Y. Gao, Phys. Lett. B655, 228 (2007), 0707.3648. 
A. Preitas, P. Schwaller, and D. Wyler (2008), 0806.3674. 

A. W. Strong et al., Astron. Astrophys. 444, 495 (2005), astro-ph/0509290. 
M. Aguilar ct al. (AMS-01), Phys. Lett. B646, 145 (2007), astro-ph/0703154. 
R. Bernabei et al. (DAMA), Eur. Phys. J. C56, 333 (2008), 0804.2741. 
O. Adrian! et al. (2008), 0810.4995. 



29 



